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Ionizing wave via high-power HF acceleration 
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Abstract. Recent ionospheric modification experiments with the 3.6 
MW transmitter at the High Frequency Active Auroral Research Program 
(HAARP) facility in Alaska led to discovery of artificial ionization descending 
from the nominal interaction altitude in the background F-region ionosphere 
by ~60 km. This paper presents a physical model of an ionizing wavefront 
created by suprathermal electrons accelerated by the HF-excited plasma 
turbulence. 



Introduction 

High-power HF radio waves can excite electrostatic 
waves in the ionosphere near altitudes where the in- 
jected wave frequency fo matches either the local plasma 
frequency f p ~ 9y / 7i^ kHz (the density n e in cm -3 ) or 

the upper hybrid resonance f uhr = yj~ff+~f% (/ c is 

the electron cyclotron frequency) [e.g., Gurevich, 1978]. 
The generated waves increase the bulk electron tem- 
perature to T e =0.3-0.4 eV, while some electrons are 
accelerated to suprathermal energies e — ^mv 2 up to 
a few dozen eV [e.g., Carlson et al., 1982; Rietveld et 
al, 2003]. Upon impact with neutrals (N 2 , 2 , O), 
suprathermal electrons excite optical emissions termed 
Artificial Aurora (AA) [e.g., Bernhardt et al, 1989; 
Gustavsson and Eliasson, 2008]. 

Heating-induced plasma density modifications are 
usually described in terms of chemical and transport 
processes [e.g., Bernhardt et al., 1989; Djuth et al., 
1994; Dhillon and Robinson, 2005; Ashrafi et al, 2006]. 
However, the Pedersen et al. [2009; 2010] discovery 
of rapidly descending plasma layers seems to point to 
additional mechanisms. Pedersen et al. [2010, here- 
after P10] suggested that the artificial plasma is able 
to sustain interaction with the transmitted HF beam 
and that the interaction region propagates (downward) 
as an ionizing wavefront. In this paper, the formation 
of such downward- propagating ionizing front is ascribed 
to suprathermal electrons accelerated by the HF-excited 
plasma turbulence. 



Ionizing wave 

The descending feature is evident in Figure [TJ which 
is representative of P10 Fig. 3 with the regions of ion- 
line (IL) radar echoes from the MUIR incoherent scatter 
radar (courtesy of Chris Fallen) overlaid [c.f. Oyama et 
al., 2006]. Shown are sequential altitude profiles of the 
green-line emissions (A =557.7 nm, excitation poten- 
tial e g «4.2 eV) observed at the HAARP facility on 17 
March 2009. Here, the O-mode radio beam was injected 
into the magnetic zenith (MZ), i.e. along the magnetic 
field Bo, at the effective radiative power Po[MW] «440 
and frequency f =2.85 MHz (2/ c at h 2fc «230 km). 



The contours of f p = fo or n e 
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and f u hr = fo (violet) are inferred from ionograms ac- 
quired at 1 min intervals [P10 Fig. 3]. The regions 
of enhanced IL are shown in green color. Note, the 
blue-line emissions at 427.8 nm (not shown) coincided 
with the green-line emissions, as seen looking from the 
HAARP site [P10]. 

During the first 2 min in the heating, the artificial 
plasma is confined to the bottomside of the F layer at 
altitudes h >180 km. The corresponding descent of the 
IL scatter is similar to that described by Dhillon and 
Robinson [2005] and Ashrafi et al. [2006]. A sudden 
brightening of AA and increased speed of descent of the 
artificial plasma 'layer' (patch) in the HF-beam center 
occurs near 180 km, while its peak plasma frequency 
foFa reaches fo. In fact, the optical data shows [P10] 
that this patch is fairly uniform near ^180 km but then 
becomes a ^20- km collection of (||Bo) filaments a few 
km in diameter. While the degree of inhomogeneity 
of the descending patch increases, its speed, V b s —0.3 
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Figure 1. Time-vs-altitude plot of 557.7 nm optical 
emissions (black color) along Bo- Blue (violet) lines 
indicate the matching altitudes /o = f p (f u h)- The 
dashed line indicates ft-2/ c - The transmitter on periods 
are indicated. Shown in green is the MUIR IL intensity 
(courtesy of Chris Fallen). Horizontal blips are stars 
passing through the view. 

km/s, appears to be constant until «160 km. Then, the 
artificial plasma slows down, staying near the terminal 
altitude h m i n wl50 km before the emissions retreat in 
altitude near the end of 4-min injection pulse. During 
a continuous 'on' period, the artificial plasma near /i m i n 
was quenched several times, initiating the process over 
again from higher altitudes. 

Hereafter, we focus on the descending feature at 
h <180 km where foFa > fo or n e > n c . Enhanced 
427.8-nm emissions indicate the presence of electrons 
with energies e > £b ~18.7 eV, exceeding the ionization 
energies of iVjj, O2, and O. The ionization rate q a is 
given by 

dn e /dt = q a = n a ■ (i/ ion (e)) (1) 

where Vi 0n is the ionization frequency, and (...) means 
averaging over the accelerated distribution of the den- 
sity n a . Hereafter, we employ Majeed and Strickland's 
[1997] electron impact cross-sections and the Hedin 
[1991] MSIS90 model for the densities [N 2 ], [O], and 
[0 2 ] on 17 March 2009. 

At each time step ti, artificial ionization occurs near 
the critical altitude h c (ti), defined from the condition 
n e (h c ) = n c . The density profile just below h c is repre- 
sented as follows 



n e (x,ti) — n c ■ \T/ (x) 



(2) 



Here x = £/Ln, £ = (h c — h) / cosao is the distance 
along Bq, ao is the conjugate of the magnetic dip angle 



(«15° at HAARP), and L\\ is the (||B ) extent of the 
ionization region. 9(x) is a monotonic function satis- 
fying the conditions \I/ (0) >1 and <Cl at x >1, 
since the ambient plasma density no <C n c at h <180 
km. As the ratio 6 e (e) of inelastic (vu) to elastic (v e i) 
collision frequencies is small, the accelerated electrons 
undergo fast isotropization due to elastic scattering and 

thus L|| ~ (liony/Se/2^, where l ion = v/v lon [c.f. Gure- 
vich et al, 1985]. 

Evidently, as soon as at some point Xi < 1 the density 
n e (xi, ti+At) ~ q a (£i)-At reaches n c , the critical height 
shifts to this point, i.e. h c (ti+\) ~ h c (ti) — Lu -x. These 
conditions define the ionization time, ~ q a /n c , and 
the speed of descent 



V d = \dh c /dt\ 



L \\ T ion 



vy/ 5e/2)n a /n c (3) 



Note, eq. ([3]) contains no dependence on the total 
neutral density N n and hence predicts Vd ~ const(/i), if 
the same distribution of accelerated electrons is created 
at each step. As (^S^ 2 v^ ~1.5T0 6 m/s, we get from 

cq. <j3j) that the value of Vd (|3]) matches V bs at n a = 

60 cm" 3 . 
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Discussion and conclusions 

We now turn to justify this acceleration-ionization- 
descent scenario. Enhanced IL echoes, like in Fig- 
ure [1] usually result from the parametric decay in- 
stability (PDI;) and oscillating two stream instability 
(OTSI) of the pump wave near the plasma resonance 
[e.g., Mj0lhus et at, 2003]. The latter develops if the 
relative pump wave energy density Wo = | Eq | / %im c T e 

exceeds Wth — — \~ -j^p, where k is the plasma wave 
number, vj- is the collision frequency of thermal elec- 
trons, and L^ 1 = |Vlnn e |. The free space field of the 
pump wave is Ef s ~5.5\/Po/r w0.65 V/matr =180 km 
(at the HF-beam center) or W fs ~5-10" 4 at T e =0.2 eV. 
For incidence angles 6 < arcsin (^/ f c / (f c + f ) sinao^j , 
the amplitude in the first Airy maximum is Ea ~ 
(2tt/ sinao ) 2/3 (/oin/c) 1/6 -B/ s [e.g., Mj0lhus et al, 2003] 
or Wa ~ 0.1(L„/L ) 1 / 3 , where L —30 km. For injec- 
tions at MZ, following Mj0lhus et al. [2003] one gets 
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As W\ [i (the electron-to-ion mass ratio), 



we get k 



1/4 



[e.g., Alterkop et al. 



1973] or krjj ~ 1/40 (rjj is the Debye radius) yield- 
ing Wth — 10~ 4 (io/£„). The 'instant' gradient-scale 
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of the artificial layer is L n ~3— s-1 km at 180— > 150 km 
(see below) gives W th «(1^3)-10- 3 . Thus, OTSI can 
easily develop in the first Airy maximum. In turn, PDI/ 
can develop in as many as ~30 Airy maxima over a dis- 
tance l a ~1 km [c.f. Djuth, 1984; Newman et ai, 1998]. 
At T e /Ti < 4, PDI; is saturated via induced scatter- 
ing of Langmuir (I) waves, piling them up into 'wave 
condensate' (k — > 0) [e.g., Zakharov et ai, 1976]. The 
condensate is subject to OTSI, thereby leading to strong 
(cavitating) turbulence and electron acceleration [e.g., 
Galeev et ai, 1977]. 

At Wi/noT e < (fc/fp) 2 , the acceleration results in a 
power-law (||Bo) distribution ate max > en = ^mu 2 > 
£min [Galeev et ai, 1983; Wang et al., 1997] 

Fj(£||) ~ n a {2p a - l)/v min ■ (e mi „/£||) Pa (4) 

where p a ~0.75-l. The density n a and e m i n are deter- 
mined by the wave energy Wi trapped by cavitons and 
the joining condition with the ambient electron distri- 
bution Fi'(e min ) = F (e m in). If F is a Maxwellian dis- 
tribution, this gives £™ in ps 10T e and n™ s» 10~ 4 n e . 
When background suprathermal (s) electrons of the 
density n s are present, then F (s >• T e ) — > F s (e) and 
e min ~ 30 {n s T e /Wif /5 eV [e.g., Mishin et ai, 2004], 
yielding e m j n < 10 eV at n e = n c , Wi ~ 10 -3 , and 
n s <10 cm -3 . In the ionizing wave, a natural source of 
the s-electrons is ionization by those accelerated elec- 
trons that can propagate from £ ~ to £ ~ Ln (see 
Figure HJ. 

We can now evaluate the excitation and ionization 
rates. The column 427.8-nm intensity in Rayleighs (R) 
is given by 

/ « W~ 6 A b JdtJ a b (e)^ a {e,0de ■ [N 2 (0] (5) 

Here a b is the excitation cross section of the N^i^N) 
state, A b «0.19, $ a = -=%-F a is the differential num- 
ber flux, and F a (e) ~ na ER ^§^v~f n {e min /e) Pa+1 is an 
isotropic distribution to which the accelerated distribu- 
tion Fa 10} is transformed at distances |£| > v/v e i due 
to elastic scattering [c.f. Gurevich et ai, 1985]. 

Integrating eq. §5§ over the energy range e b < s <10 2 
eV at p a =0.85 yields the brightness of a A/i-km column 
AI(h c ) « 2.5 • 10~ 12 : n a [N 2 {h c )] ■ Ah R near altitude h c , 
given that Ah <C H n ~8 km (the atmosphere scale- 
height). The total intensity / is defined by the vertical 
extent of the (excitation) layer A b , where e(£) > e b . It 
can be evaluated using the Majeed and Strickland [1997] 
loss function L(e) = J2j^j( £ ) with j designating A?2, 
2 , and O. Outside the acceleration layer, i.e. |£| > l a , 



altitude, km 

2 150 200 250 




20 40 60 80 100 

energy, eV 



Figure 2. (a) Altitude profiles e(e ,£.) at e =10, 
15, ...100 eV and ho =150, ...200 km. (b) Half-widths 
A g (thin lines) and Ab (thick) of the green- and blue- 
line excitation layers near h c =160 (circles) and 180 
(solid lines) km. 

the energy of an electron of the initial energy £0 at a 
distance £ from the origination point ho is 

rh +i 

£(£ ,0-£o-/ L(e(z)W2/S e (e(z))dz (6) 

Jho 

Figure [5^ presents the results of calculations of eq. 
© for £ =10, 15, ... 10 2 eV and h =150, 160, 180, 
and 200 km. The altitude profiles at e >5 eV and hence 
the layers of excitation/ionization are nearly symmet- 
ric about hg at h <180 km. Panel b shows the half- 
widths A g and Ah of the green- and blue-line excita- 
tion layers about h c =180 and 160 km as function of 
£0. The half- width of the ionization layer Ai on (not 
shown) is s» Ah at £0 >20 eV. Since A b < H n , we can 
estimate the 427.8-nm intensity at h c =180— >160 km 
as I\ he ~5-10" 12 n a • [N 2 {h c )] (A b (h c )) ~(0.16^0.2)-n Q 
R. Comparing JL with the spatially-averaged intensi- 
ties % sal0^5 R [P10 Fig. 1] yields n a ~60^25 cm- 3 . 
Note that n a ~ n a at 180 km, in agreement with a uni- 
form structure, while spatial averaging underestimates 
n a inside the ^km-scale filaments at 160 km. 

Calculating the ionization frequency in eq. (QJ with 
F a {e) gives {p ion ) w n* on ■ {[N 2 ] + \[0]+ 0.95[O 2 ]) s"\ 
where n* on = (v<Ji on ) /n a «1.8-10~ 8 cm 3 s~ 1 is the 
coefficient of ionization of N 2 . The total ionization 
rate qa ~ 10 4 cm~ 3 s~ 1 greatly exceeds recombination 
losses s» 10~ 7 ri 2 ps 10 3 cm _3 s~ 1 (the main ion compo- 
nent at these altitudes is NO + ). This justifies the use 
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Figure 3. The ionization coefficients of N2, O2, and 
and the excitation coefficient of the A^ l "( 1 iV) state vs. 

^max- 

of eq. (CL} for evaluating the artificial plasma density. 
Taking an average energy loss per ionization ~20 eV 
results in the column dissipation rate <0.1 mW/m 2 or 
<10% of the 440- MW Poynting flux, consistent with 
PlO's estimates. 

As shows Figure 02 the coefficients of ionization and 
blue-line excitation by accelerated electrons decrease by 
a factor of -2 (10) between e max =10 2 and 50 (30) eV. 
The Liouville theorem predicts F(e — Ae(£o, £), ho + 
= Fo(so, h ), where Ae(e ,C) is given by the integral 
in eq. ([6]). Thus, the gradient scale-length L n of the 
artificial plasma is about the distance £ 50 , defined by 
the condition Ae(10 2 ,£ 50 ) ~50 eV. Numerically, we get 
£ 50 ~ A 6 (50) or L n «3->d.5 km near h c =180^160 km 
and q a L n /n c ~ V bs> as predicted by eq. Note that 
the artificial plasma density profiles derived from iono- 
grams indeed have ~l-km gradient scale-lengths near 
150 km [c.f. P10 Fig. 2]. 

Figure [Q shows that the descent slows down below 
160 km and ultimately stops at ft. m i n «150 km. The 
presence of IL and bright green-line emissions indicate 
that plasma turbulence is still excited and efficiently ac- 
celerates electrons above 4 eV. However, the blue-line 
emissions almost vanish [P10], thereby indicating only 
few accelerated electrons at e > £;,. That this is in no 
way contradictory follows from the fact that inelastic 
losses increase tenfold between 10 and 20 eV. Acceler- 
ation stops at e = £ max <§Cl00 eV when vu{e max ) ex- 
ceeds the acceleration rate mD||(M max )/87re max , where 
J 1 1 2 

D || (u) ~ in v mu \EkA and fey = lo p /u [Volokitin and 
Mishin, 1979]. The critical neutral density is roughly 
estimated as ~5T0 n cm -3 , i.e N n at ~150 km. The 
fact that the artificial plasma stays near h m - ln indi- 
cates that ionization is balanced by recombination or 
<7™ m ~ 10~ 7 n 2 w 0.lqi d \ which at n a ~ corre- 



sponds to e max sa30 eV (Figure [3]). 

A mechanism for generating km-sized filaments be- 
low 180 km could be the thermal self-focusing insta- 
bility (SFI) near h c , resulting in a broad spectrum of 
plasma irregularity scale sizes [e.g., Guzdar et ai, 1998]. 
Significantly, ~km-scale plasma irregularities grow ini- 
tially but within 10s of seconds thermal self-focusing 
leads to smaller (10s to 100s meters) scale sizes. Dur- 
ing descent, the critical altitude moves downward by 
several km within 10 s, thereby precluding further de- 
velopment of SFI, while the ~km-scale irregularities 
have sufficient time to develop. When the descent rate 
drops, small-scale irregularities can fully develop and 
scatter the HF beam, thereby impeding the develop- 
ment of OTSI/PDIi and hence ionization. As soon as 
the artificial plasma decays, SFI falls away and hence 
irregularities gradually disappear. Then, the artificial 
plasma can be created again. This explains why the 
artificial layer ceases and then reappears (Figure [1} . 

In conclusion, we have shown that the artificial 
plasma sustaining interaction with the transmitted HF 
beam can be created via enhanced ionization by suprather- 
mal electrons accelerated by Langmuir turbulence near 
the critical altitude. As soon as the interaction region 
is ionized, it shifts toward the upward-propagating HF 
beam, thereby creating an ionizing wavefront, which 
resembles Pedersen et al's [2010] descending artificial 
ionospheric layers. 

Acknowledgments. 

This research was supported by Air Force Office of Sci- 
entific Research. We thank Chris Fallen for providing the 
MUIR IL data. 

References 

Alterkop, B., A. Volokitin, V. Shapiro, and V. 
Shevchenko (1973), Contribution to the nonlinear 
theory of the "modified" decay instability, J ETP Let- 
ters, 18, 24. 

Ashrafi, M., M. Kosch and F. Honary (2006), Heater- 
induced altitude descent of the EISCAT UHF ion- 
line enhancements: Observations and modeling, Adv. 
Space Res., 38, 2645. 

Bernhardt, P., C. Tepley, and L. Duncan (1989), Air- 
glow enhancements associated with plasma cavities 
formed during ionospheric heating experiments, J. 
Geophys. Res., 94, 9071. 

Carlson, H., V. Wickwar, and G. Mantas (1982), Ob- 
servations of fluxes of suprathermal electrons accel- 



MISHIN AND PEDERSEN: HF-INDUCED IONIZATION 



5 



erated by HF excited Langmuir instabilities, J. Atm. 
Terr. Phys., 12, 1089. 

Djuth, F., (1984), HF-enhanced plasma lines in the 
lower ionosphere, Radio Sci., 19, 383. 

Djuth., F., P. Stubbe, M. Sulzcr, H. Kohl, M. Ri- 
etveld, and J. Elder (1994), Altitude characteris- 
tics of plasma turbulence excited with Troms0 su- 
perheater, J. Geophys. Res., 99, 333. 

Dhillon, R. S., and T. R. Robinson (2005), Observa- 
tions of time dependence and aspect sensitivity of 
regions of enhanced UHF backscatter associated with 
RF heating, Ann. Geophys., 23, 75. 

Galeev, A., R. Sagdeev, V. Shapiro, and V. Shevchcnko 
(1977), Langmuir turbulence and dissipation of high- 
frequency energy, Sov. Phys. JETP, 46, 711. 

Galeev, A., R. Sagdeev, V. Shapiro, and V. Shevehenko 
(1983), Beam plasma discharge and suprathermal 
electron tails, in Active Experiments in Space (Alp- 
bach, Austria), SP-195, pp. 151, ESA, Paris. 

Gurevich, A., Y. Dimant, G. Milikh, and V. Vaskov 
(1985), Multiple acceleration of electrons in the re- 
gions high-power radio-wave reflection in the iono- 
sphere, J. Atmos. Terr. Phys., ^7, 1057. 

Gustavsson, B., and B. Eliasson (2008), HF radio wave 
acceleration of ionospheric electrons: Analysis of 
HF-induced optical enhancements, J. Geophys. Res., 
113, A08319, doi:10.1029/2007JA012913. 

Guzdar, P. N., P. K. Chaturvedi, K. Papadopoulos, and 
S. L. Ossakow (1998), The thermal self-focussing in- 
stability near the critical surface in the high-latitude 
ionosphere, J. Geophys. Res., 103, 2231. 

Hedin, A. (1991), Extension of the MSIS thermospheric 
model into the middle and lower atmosphere, J. Geo- 
phys. Res., 96, 1159. 

Majeed, T., and D. J. Strickland (1997), New survey of 
electron impact cross sections for photoelectron and 
auroral electron energy loss calculations, J. Phys. 
Chem. Ref. Data, 26, 335. 

Mishin, E., W. Burke, and T. Pedersen (2004), 
On the onset of HF-induced airglow at mag- 
netic zenith, J. Geophys. Res., 109, A02305, doi: 
10.1029/2003JA010205. 

Mj0lhus, E. Helmersen, and D. DuBois (2003), Geomet- 
ric aspects of HF driven Langmuir turbulence in the 
ionosphere, Nonl. Proc. Geophys., 10, 151. 

Newman, D., M. Goldman, F. Djuth, and P. Bern- 
hardt (1998), Langmuir turbulence associated with 
ionospheric modification: Challenges associated with 
recent observations during a sporadic-E event, in: 
Phys. of Space Plasmas., ed. by T. Chang and J. 
Jaasperse, v. 15, p. 259, MIT, Cambridge, MA. 



Oyama, S., B. J. Watkins, F. T. Djuth, M. J. 
Kosch, P. A. Bernhardt, and C. J. Heinselman 
(2006), Persistent enhancement of the HF pump- 
induced plasma line measured with a UHF diagnos- 
tic radar at HAARP, J. Geophys. Res., Ill, A06309, 
doi:10.1029/2005JA011363. 

Pedersen, T., B. Gustavsson, E. Mishin, E. MacKenzie, 
H. C. Carlson, M. Starks, and T. Mills (2009), Opti- 
cal ring formation and ionization production in high- 
power HF heating experiments at HAARP, Geophys. 
Res. Lett, 36, L18107, doi:10.1029/2009GL040047. 

Pedersen, T., B. Gustavsson, E. Mishin, E. Kendall, 
T. Mills, H. C. Carlson, and A. L. Snyder (2010), 
Creation of artificial ionospheric layers using high- 
power HF waves, Geophys. Res. Lett., 37, L02106, 
doi:10.1029/ 2009GL041895. 

Rietveld, M., M. Kosch, N. Blagoveshchenskaya, V. Ko- 
rnienko, T. Leyser, and T. Yeoman (2003), Iono- 
spheric electron heating, optical emissions and stria- 
tions induced by powerful HF radio waves at high lat- 
itudes: Aspect angle dependence, J. Geophys. Res., 
108, 1141, doi: 10.1029/2002JA009543. 

Wang, J., D. Newman, and M. Goldman (1997), Vlasov 
simulations of electron heating by Langmuir tur- 
bulence near the critical altitude in the radiation- 
modified ionosphere, J. A. S. -T. P., 59, 2461. 

Volokitin, A., and E. Mishin (1979), Relaxation of an 
electron beam in a plasma with infrequent collisions, 
Sov. J. Plasma Phys. 5, 654. 

Zakharov, V., S. Musher, and A. Rubenchik (1976), 
Weak Langmuir turbulence of an isothermal plasma, 
Sov. Phys. JETP, 42, 80. 

Space Vehicles Directorate, Air Force Re- 
search Laboratory, Hanscom AFB, MA 
01731 (e-mail: Evgeny.Mishin@hanscom.af.mil; 
Todd . Pcdersen@hanscom . af . mi) 



This preprint was prepared with AGU's WTpft. macros v4, 
with the extension package 'AGU++' by P. W. Daly, version 1.6b 
from 1999/08/19. 



